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ANALYSIS OF TEMPERATURE DISTRIBUTION IN LIQUID-COOLED TURBINE BLADES ' 


By Joux N. B. Livincoop and W. Brron Brown 


SUMMARY 


The temperature distribution in liquid-cooled turbine blades 
determines the amount of cooling required to reduce the blade 
temperature to permissible ralues at specified locations. This 
report presents analytical methuds for computing temperature 
distributions in liquid-cooled turbine blades, or in simplified 
shapes used to approximate sections of the blades. 

The individual analyses are first presented in terms of their 
mathematical development. By means of numerical examples, 
comparisons are made between simplified and more complete 
solutions and the effects of several variables are examined. 
Nondimensional charts. to simplify some temperature- 
distribution calculations are also given. 

It was found that for blade materials having low thermal 
conductimties (corresponding to blade materials in current use) 
a one-dimensional spanwise temperature-distribution equation 
ig applicable near a coolant passage. The cooled part of the 
blade, irrespective of blade length, was found to remain at a 
determinable uniform temperature, called the prevalent blade 
temperature. The prevalent blade temperature increased abvuut 
150° F for an increase in the difference between the effective gas 
temperature and the coolant temperature of 1000° F for a range 
of effective gas temperature from 2000° to 5000° F. 

In cases where rim cooling is insignificant, it was found that a 
one-dimensional chordwise temperature-distribution equation 
was eufficiently near the relaxation solution for the actual blade 
shape for a first-order approximation. 


INTRODUCTION 


A limitation on design and performance of aircraft gas- 
turbine power plants is the strength of the turbine materials, 
which decreases as temperature increases. This limitation 
may be greatly alleviated, even when nonstrategic materials 
are used, by the application of turbine cooling. An extensive 
study of both liquid and air cooling of turbine blades has been 
carried out at the NACA Lewis laboratory since 1945. 

The temperature distribution in a turbine blade determines 
the amount of cooling required to reduce the blade tempera- 
ture to permissible values at specified locations. Analytical 
studies have been conducted to predict temperature distribu- 
tions in cooled turbine blades. A summary of the analytical 
methods developed for calculating spanwise temperature 
distributions for three types of air-cooled turbine blade is 
presented in reference 1. The present report summarizes 
the analytical methods developed for calculating temperature 
distributions in liquid-cooled turbine blades. Parts of these 
analyses are presented in references 2 and 3. 


Because the calculation of a generalized three-dimensional 
temperature distribution through a liquid-cooled turbine 
blade is too tedious to be of practical value and because some 
knowledge of the distribution in a radial direction is neces- 
sary, the first investigation presented herein is limited to a 
one-dimensional spanwise distribution. 
however, is valid only in the immediate vicinity of the 
coolant passages. In order to study more accurately the 
temperature distribution in the trailing section of the blade, 
a three-dimensional solution is determined for a rectangular 
parallelepiped used to approximate a blade trailing section. 
This solution serves to determine the importance of blade 
thermal conductivity and of distance from the coolant. pas- 
sage. Because numerical calculations for a rectangular 
parallelepiped produce & constant spanwise temperature in 
the region of the blade beyond the influence of rim cooling 


This distribution, — 


-— — — u oe 


(that is, in the region of the blade unaffected by conduction 


to the rim), one-dimensional chordwise temperature distribu- 
tions are then calculated for shapes approximating various 
parts of a blade cross section; these investigations are valid 
for the part of the blade beyond the influence of rim cooling. 
Finally, a two-dimensional temperature distribution through 
the cross section of an actual blade shape is obtained in order 
to determine the accuracy of the previously determined 
approximate solutions and to study the effect on temperature 
distribution of a peripheral variation of gas-to-blade heat- 
transfer coefficient. 
Numerical examples based on coefficients available from 
unclassified sources are included. 


METHODS OF ANALYSIS 


One-dimensional spanwise, three-dimensional, one-dimen- 
sional chordwise, and two-dimensional chordwise temperature- 
distribution equations for liquid-cooled turbine blades 
or for simplified shapes used to approximate sections of 


ee oe 


liquid-cooled turbine blades are presented. For simplicity 
in presentation, the individual analyses are first formulated 
in terms of the mathematical methods. The comparative 
applicabilities of the methods of analysis are established in 
later sections of the report. 

For all the analyses presented, the following conditions are 
assumed: 

(1) The coolant temperature is constant at the average 
temperature of the coolant, or the coolant forms a constant- 
temperature reservoir. This assumption is valid when the 
change in coolant temperature is small in comparison with 


the temperature difference between the gas and the coolant. 


! Supersedes NACA TN 2821, “Analysis of Temperature Distribution in Liquid-Cooled Turbine Blades” by John N. B. Livingood and W. Byron Brown, 1951. 
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(2) The blade-to-coolant heat-transfer coefficient is con- 
stant. 

(3) The thermal conductivity of the blade is constant. 

(4) The effective gas temperature at all parts of the blade 
is the same. 

(5) Radiation effects are considered to be included in. the 
heat-transfer coefficients. The following individual analyses 
are presented: 

One-dimensional spanwise temperature distribution —A 
one-dimensional spanwise calculation gives a radial tempera- 
ture distribution valid in a part of the blade near a coolant 
passage. In this particular case, the radial temperature 
distribution is carried through the blade and through the 
turbine rotor. Heat transfer between the turbine rotor and 
the fluid on the outside of the rotor is assumed to take place 
and a constant value of gas-to-blade heat-transfer coefficient 
is assumed. For simplicity, a blade of constant cross- 
sectional area and perimeter is considered. 

For blades with little taper, the method is applicable if 
average cross-sectional area and perimeter are considered. 
In general, however, a numerical solution is necessary for 
tapered blades. 

Three-dimensional temperature distribution through sim- 
plified trailing section.—A. three-dimensional temperature 
distribution for the trailing section of the blade is obtained 
by considering this part of the blade as a rectangular 
parallelepiped. The blade-root temperature 1s considered 
constant, the gas-to-blade heat-transfer coefficient is constant, 
and it is assumed that no heat flows across the median plane 
of the simplified blade section. 

An approximate solution for tapered blades is obtainable 
by this method if the blade is considered in small sections, 
average values of cross-sectional area and perimeter are 
considered for each section, and appropriate compatibility 
relations are satisfied. 

One-dimensional chordwise temperature distribution 
through simplified shapes.—Rectangular and trapezoidal 
sections, concentric circle annuli, and sections between 
parallel plates are considered as the simplified shapes used 
to approximate parts of a turbine-blade cross section. 
The gas-to-blade heat-transfer coefficient is. considered 
constant. For comparison purposes, the rectangular and 
trapezoidal sections are constructed so that lengths and 
areas are identical. 

Two-dimensional chordwise temperature distribution.— 
Numerical two-dimensional temperature distributions 
through the cross section of an actual blade shape are de- 
termined. Constant and variable gas-to-blade heat-transfer 
coefhicients are considered. 


ONE-DIMENSIONAL SPANWISE TEMPERATURE DISTRIBUTION 


The spanwise temperature-distribution equations for a 
liquid-cooled turbine, a section of which is shown in figure 1, 
are derived in reference 2 and are reproduced in appen- 
dix B. In this spanwise case, the investigation was extended 
from blade tip to rotor hub. 

For convenience, the turbine was divided into four 
sections (fig. 1), and the temperature-distribution equations 
were obtained from heat balances for differential elements 
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__ FIGURE 1.—Arrangement of internel-cooling pessages In turbine. 


in each section. Inasmuch as chordwise conduction was 
neglected, the following equations, valid only in the neighbor- 
hood of the coolant passages, were obtained. For the 
uncooled section of the blade (fig. 1, section 1), it was found 
that 


Ts, Tg e — Ci cosh A (i+ Ca) (1) 
where C, and C; are integration constants and 
= helo rt) 
kpAp 1 


(AH symbols are defined in appendix A.) The value of C 
is determined by use of the boundary condition at the blade 
tip (where heat enters by convection only). 





—kpAp 1 (52) =, Ag (T; e— Ts, r) 
1/T 
or 
: hotp 1 $ Meh sie: fe 
tanh AC peh =e) 


For the cooled section of the blade (fig. 1, section 2), it 
was found that 


Tg = Tege p— Cet — Oei (2) 
where C, and (, are integration constants, 


holo Ty, e-t As, alt, ol Ps 


Ta o= daha ee 
and 
hely-t he h N$ 
wa eae 
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The prevalent blade temperature Ts, is the temperature 
the blade would assume if no heat flows in or out of the blade 
ends. 

For the rim section (fig. 1, section 3), with assumed con- 
stant area, 


Tn = G+ Cse Cee- Er (3) 


where C, and Cs are integration constants, 


G= 





Ears ashta. elta T +(hi, alt. mal 
Aarau e+ (h, att, aa 


E= tars ae tets ali, weasel 


Aim ch rim 


and J, is the fluid temperature on the outside of the rotor 
and (her alss)ae is the average of the values of hiala found 
for the blade coolant passages, the rim circumferential 
passages, and the radial passages through the rotor (fig. 1). 

For the rotor section (fig. 1, section 4), with assumed con- 
stant strength, 


=J- C; cosh H, (4) 


where C; is an integration constant, and where 


dor, alela +(h, ale, ‘aol 


—— eae | 


J= vt, 4 
AIMS as hat (hi, ilia sav 





and 


H= rr, ar hat (hi, als, let Rede | 


Kelp 


Values of the six intergration constants Ci, Cz, Cy. Cs, Co, 
and Cy are found by solving simultaneously the six equations 
resulting from equating temperatures and heat flows at the 
various Junction points of the four sections of the turbine. 
(Subscripts a and 8 used after numerical subscripts designate 
the end nearest the blade tip and the end nearest the rotor 
hub of the sections to which the numerical subscripts refer, 
respectively.) At the junctions of sections 1 and 2 


1 C: cosh A(zıs+ C) = Ta, p— Cye 


(5) 
ke sip, 1pC1A sinh A(x Lat C)=BkpAp > (Ces man 
Ceara) hi (Ap 1 p— Az, > a) (To, p— Oe >a 


( e hee 7) (6) 


Equation (6) equates the heat leaving section 1 to the sum 


of the heat entering the metal of section 2 and the heat 
entering the coolant at that part of the blade where the 
inlet and outlet passages are connected. An approximation 
in the procedure has been introduced at this point because 
a separate blade section for the part of the blade containing 
the connecting passage between inlet and outlet passages is 
not introduced. A numerical calculation showed that use 
of such a section would slightly decrease the temperature at 
the blade tip. 
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At the junction of sections 2 and 3, 


T's, p Cert. — Oe 72.6 — G + Cresa L Cie” ™s. (7) 
kgn, »B (Cren S Ce 72. 6) + 
Aue Leas aon Tp pt Cae 284 
“Lp, 
Ce" **2.8)— Ekg Ag 2 Fe (Cee es a— Ce" ý a) (8) 


Equation (8) equates the sum of the heat leaving section 
2 and that entering section. 3 directly from the hot gases to 
the total heat entering section 3. 

At the junction of sections 3 and 4, 


G+ Cye™30+C,e7 e40; cosh Hra 


kpE(—Cye ast Oye nl) 2 ksHC; sinh Hra 





(10) 


THREE-DIMENSIONAL TEMPERATURE DISTRIBUTION THROUGH 
SIMPLIFIED TRAILING SECTION 

The spanwise temperature-distribution equations previ- 
ously presented are valid only in the neighborhood of the 
coolant passages. Most currently used turbine blades are 
so shaped that coolant passages cannot be located near the 
trailing edge. Because of the impossibility of placing reason- 
ably large coolant passages near the trailing edges of con- 
ventional turbine blades and because the trailing sections 
seem most likely to be the hottest parts of the blade, other 
detailed studies were made of temperature distributions. 
The first study was devoted to approximating the trailing 
section by a rectangular parallelepiped and to determining 


the three-dimensional temperature distribution through this 


parallelepiped (fig. 2). 
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FiGUBRE 2.— Rectangular parallelepiped used to approximate trailing section of 4-Inch turbine 
blade for three-~imensfonal analysis. 


(9) 
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The median plane of the rectangular parallelepiped was 
chosen as the plane 2=0, and symmetry was assumed about 
this plane (no heat flow was assumed to cross this plane, 
as previously stated). Boundary conditions at the blade 
trailing edge and blade tip were simplified by assuming the 
blade width and blade length to be extended by a distance 
equal to one half the blade thickness. In figure 3, the edge 
MN at temperature T gained some heat from the gas stream. 
The extended surfaces NO and MP were at nearly the same 
temperature T and no heat entered the edge OP; these 
surfaces therefore gamed practically the same amount of 
heat as the actual exposed edge. The validity of this ap- 


proximation is discussed in reference 4 (pp. 216 and 217). 
Distances increased by 7/2 are denoted by primes. 

The derivation of the three-dimensional temperature dis- 
tribution, originally derived in reference 3, is reproduced in 
appendix C. The differential equation, in final form, for 
the temperature distribution is found to be 


076 , O76 , 078 


arat dyt vat? (11) 


where 2’ and y’ denote x and y increased by 7/2, respectively. 
A solution of equation (11) satisfying the boundary condi- 
tions at x =y =z=0 Is | 


(= =>} È (Kin, cos Laz’ cosh M m, n Y’ cos Nmz -+ 


mml ne] 


Om. n cosh P, ,z’ cos Q,y’ cos Nma) (12) 


M 
Liquid Y hee where K, L, M, N, O, P, and Q are constants. Relations 
coolant te 7 among the constants are 
AOE 2 pa M.. 2=N.2-+L.? 
N 0 ma =N m T Lin anei Bipi (13) 
FiGtRE 3.—Correction for heat received by trailing edge. P m, = N a + Q n d 
Values of Km, „ and Om.» are given by the relations 
Ni. = 
x et UNM SAT, =T) ( Y 
r (2n— 1) N,2 sin 2N, 2 ky My, sinh Mm, = +h, cosh Maa) 14) 
1+ — 
oN I 
2N; 5 
and 
0, EE e E T MZ E E: p 
ma 7 = 
Qa I (11528 2@nJ +e. sin 2N, : cosh P,, nb (15) 
1+ s 
2N m 5 


Application of the boundary conditions at 2’=0’, 
N, respectively: 


tan Qi = pear : z= = eee . : -- =- — 








and 
het 
tT he 2kg | | 





It is at once obvious that from equation (16) many values 
of LE, result. Because of the periodicity of the tangent 
function, equations (17). and (18) likewise have many 
solutions. 

Values for all the constants in equation (12) can thus be 
found and the temperature can now be computed at any 
point in the rectangular parallelepiped. In reference 3, 
appendix B, it is shown that sufficiently accurate results can 
be obtained by using m=1 only. 


y =)’, and z=7/2 lead to the following methods of evaluating L, Q, and 


-e 08) 


hi? 
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ONE-DIMENSIONAL CHORDWISE TEMPERATURE DISTRIBUTIONS 
THROUGH SIMPLIFIED SHAPES 

Because the spanwise temperature distributions are valid 
only in the neighborhood of the coolant passages and because 
the three-dimensional approximate solution resulted in a 
constant spanwise blade temperature in the part of the blade 
beyond the influence of rim cooling, one-dimensional chord- 
wise temperature distributions were determined for sections 
of a liquid-cooled turbine blade that can be approximated 
by simple shapes (fig. 4). Rectangular and trapezoidal 
approximations for blade trailing sections were considered. 
Trapezoidal sections may also be used to approximate the 
leading section in some blades. In addition, analyses were 
made for leading sections approximated by concentric circle 
annuli and for the sections of blades with very little metal he- 
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FIGURE 4.—8implified shapes used in one-dimensional chordwlse analysis. 


tween blade outer surface and coolant passages approximated 
as regions between parallel plates. The temperature distri- 
butions through these simplified shapes, derived in appendix D 
and valid in regions beyond the influence of rim cooling, follow. 
Rectangular trailing section of blade removed from 
influence of rim cooling.—For a blade trailing section 
approximated by a rectangle (fig. 4(a)) the temperature- 
distribution equation is 
hy 
kp (T,,e—T;) cosh gy’ 
9 =— i l *19) 
y sinh ag cosh 97’ 


where 


é= T; s— Tg 


‘2he\t 
a- 


Trapezoidal trailing section of blade removed from in- 
fluence of rim cooling.—For a trapezoidal approximation to 
the blade trailing section (fig. 4(b)), the temperature distri- 
bution is 


and 





Aita ( 
_ 2K ks 


inair- iinit T 


Te Ty) (Mitr ACH HiHat) fae 
20 








K= (pty) 


ay T2 T1 
Ņ=tan ES 


T=H (it) Jit) tiir NHAT) 


and Ja, iJı, Ho, and iH, are Bessel functions. For a wedge- 
shaped section the temperature-distribution equation (20) 
reduces to 





E (Ts. — —T;) Jo(t t) 
EEEE (21) 
«sles était 


Concentric-circle-annulus approximation for blade leading 
section.— The equations expressing the Llade-metal tem- 
perature at the coolant passage wall Ts : at the leading edge 
of the blade Tz. and at any point in the annulus Ts are 
found to be (fig. 4(c)) 


( 2ks 
Tp. {= Tı+ En log, (By [erates (22) 
i 


fe Bo—Lp:) (23) 


2kp 
Ts, a= Trs -i Doho log, (2) 
i 


and 
(D) 
log 7 
Ts= Ts, H Nae (24) 
log. (5) 
where 
TaT np r 8) 
a Dy (on thi 
0 é D, 


Section of blade approximated by parallel plates.—The 
equations expressing the blade-metal temperatures at the 
coolant-passage wall (fig. 4 (d)) Ts; and at the blade outer 
edge Tg o are 


r Gi) eTa ethT) 


Hi 5) Goth) ma 
and 
_ Tet aii, =p) (hel eet hs T) - 


nEs 5) (he +h;) 


where o is the distance between the parallel plates. 
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TWO-DIMENSIONAL CHORDWISE TEMPERATURE DISTRIBUTION 


The temperature at any point on the blade cross section 
may be found by solving the Laplace differential equation 
given in terms of the temperature difference 0 


O76 , O76 
ator? 


where X and Y are the Cartesian coordinates in the blade 
cross section. The boundary conditions, expressed in 
terms of the partial derivative of 0 in the direction normal 
to the boundary, are 


(28) 


ob ho 
on ka 120) 
at the outer boundary and 
Ob hy dip 
on ke (dri T; 8) (3 0) 


at the coolant-passage boundary. 


A closed solution to equation (28) cannot be obtained 


because of the impossibility of applying the given boundary 
conditions along the odd-shaped boundary of a turbine 
blade (fig. 5). A numerical solution is available, however, 
by application of the relaxation method (reference 5). A 
sketch of the blade cross section is covered by a network of 
points (a square network was chosen). Large net spacings 
are recommended at first, and the insertion of additional 
net points may be made at any time during the solution, 
thus permitting the use of a final network of any desired 
size. Temperature estimates for the net points can be obtained 
from the solutions of the chordwise blade-temperature- 
distribution equations for approximate shapes previously 
presented and from application of the boundary conditions 
given by equations (29) and (30).. Residuals, which may 
be considered as interior heat sinks, can then be calculated 
at each net point whose immediately neighboring points 
remain within the boundary from the relation 


Qo= 61+ 89a + 5+ 0 — 4b (31) 


where the subscript zero denotes a point in the blade cross 
section and subscripts 1, 2, 3, and 4 denote the points in 
the square surrounding the point with zero subscript. 





Friaure &—Turbine-blade section showing coolant passages. 


The object of the relaxation is to reduce the values of 
Q to zero, or as close to zero as possible. When the relaxa- 
tion equation (31) is employed, the following procedure is 
used. A change in @ alters Q,, Qa, Qs, Q4 by the same change 
and Q, by minus four times this change, all other values of 6 
remaining fixed. 

Equation (31) is the finite-difference equation correspond- 
ing to the partial differential equation (28). For net points, 
some of whose immediately neighboring points lie outside the 
boundary, equation (31) must be modified. For example, if 
point 1 lies outside the boundary and point 5 is the boundary 
point between 0 and 1, the following equation applies: 


ey yaar bst +4—(3+5) fo (32) 


where d is the ratio of the distance between points 0 and 5 
to the net spacing. 

Corresponding changes in the values of Q result from the 
use of equation (32). No harm is done by overestimating 
the final values of @, as successive calculations will establish 
them again. Continued relaxation eventually reduces all 
the residuals as desired and the blade temperatures can 
finally be obtained from the definition of 8, 


é= T, en Ip 
APPLICATION OF ANALYSIS 


Typical numerical examples for a sample blade are 
presented to illustrate temperature trends and to determine 
the effects of various factors on the temperature distributions 
in a liquid-cooled turbine blade. A turbine blade whose 
external shape is similar to that of the rotor-root section of a 
conventional gas-turbine design was selected as the sample 
blade. Two 0.25-inch-diameter coolant passages were 
assumed in the blade (fig. 5), connected near the blade tip 
by a cross-over passage, and the cooling was assumed to 
occur through forced convection. No allowance was made 
in these calculations for the effects of free-convection currents 
that might be present. A gas flow of 55 pounds per second 
(equivalent to a mass velocity of about 58 lb/(sec)(sq ft)) 
and a coolant flow of 7 pounds per minute per blade were 
assumed (to insure turbulent flow and to remove the depend- 
ence of the blade-to-coolant heat-transfer coefficient on the 
length-to-diameter ratio of the coolant passage). 

For the one-dimensional spanwise case, in which the 
cooling was carried to the rotor hub, two 0.25-inch-diameter 
coolant passages running circumferentially through the rim 
and ten 0.50-inch-diameter coolant passages running radially 
through the rotor were assumed (fig. 1). For this case, 
water, ethylene glycol, and kerosene were all considered as 
possible coolants. Other calculations were made only for 
water as the coolant because of its superiority over the other 
coolants considered. 

Blade-to-coolant averago heat-transfer coefficients were 
calculated by use of formulas presented in reference 6, 
page 168 for turbulent flow and page 190 for laminar flow 
(for ethylene glycol only). Gas-to-blade average heat- 
transfer coefficients were calculated by use of formulas given 
on page 236 of reference 6. The following coefficients were 
used in the analysis: 
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Outside heat-transfer coefficient, A, 


Bin) SOIC feo s65 2 Se eee ea Seo se 222 
Btu (s60)4SG40) Ch) acount oe eines 0. 06167 
Inside heat-transfer coefficient, Ri.: 
Water, 
Biure ft) CE) essene bs cee scaseaaeice 2370 
Btu/ (sec) (sq ft) (°F)_-_------.------------ Tonine 0. 6583 
Ethylene glycol, 
Btw hr) (Sa EOC TF oc.cn cance cctee cecwe nee deceses 649 
Btu (60) (Sa OCE] nck cuccw cence coveneedere 0. 1803 
Kerosene, 
Btu (Or) (86-16) Clee es pce sees cows 510 


Btu/(see) (sq ft)(°F)---------------------------- 0. 1417 


Other numerical values used in the calculations were 
(figs. 1 and 4): 


Effective gas temperature, Taea °F_..------------------ 2000-5000 
Average coolant temperature, 7), °F_--._--.-------------- 200 
Number 0! -Dindes, Z- coche ce ccetiowee cu ecuedeneces 55 
Area, of blade cross section 1, Ai, sq ft----------- iL eee 0. 00198 
Area of blade cross section 2, Az, sq ft-------------------- 0. 0013 
Area of cross section 3, Az, sq ft_.-_...-------_-.-------- 0. 312 
Area of cross section 4, A4, 8q ft___...._._----.___.-.- _.. 0, 286 
Blade outside perimeter, le, ff...0e....-.----------------- 0. 2542 
Blade inside perimeter, las ft_..--.-.-------------.-...-. 0181 
Average radial distance, section 3, Tsen ft-......---.--.--- 0. 4917 
Average radial distance, section 4, Tiran ft_-----..-------.-- 0. 3333 
Length of chordwise trailing section, J, ft------------------ 0. 050 
Length of spanwise trailing section, b, ft.._.-...---------- 0. 3333 . 
Thickness of trapezoid at coolant passage, rz, ft._--.-.-.-.-. 0021 
Thickness of trapezoid at trailing edge, tu ft-_..--_-.------ 0. 003 
Thickness of rectangle, r, ft...-.---_-------------------- 0. 010 
Temperature at blade root, T,, °F_.....-...-.------------ 330 
Thermal conductivity, ks 
Btu (he) lt) Cae oe es ie ee ee as 15-210 
Btu/ (sec) (ft) CP) ncosooc ce en eects 0. 00417-0. 0583 
Cooling-air (in contact with rotor) temperature, Ta, °F'_..--- 0 
tleat-transfer coefficient between cooling air and rotor, he 
Bt) Hr) Gå (CF) -sedara aenn n E R 30 
Btu/ (sec) (sq MC F] -asosiari Tocana tale 0. 00874 
Average value in section 3 of l; 3h;,3 
Water, 
Btu bry (CB) 2.3 oo ct oe eee ees 30, 360 
Btu) (S60) (it) C Bosse cere i te eee 8. 433 
Ethylene glycol, 
BU ONUC] oe eee aE 8305 
Btu Gec) hO P-er a re Geet 2. 307 
Kerosene, 
Btu (MRN 2) chim een ei a na e 6533 
Bulo GOC Pekara acetate ed cees 1. 815 
Average value in section 4 of 1: sh. 4 
Water, 
Bo it) (i) COP aaa ace eeS eee cheeses 6107 
Btu/ (see) (ft) (°F) --------~---------------------- L 686 
Ethylene glycol, 
Btu? Che (66) Fonar aeaa a 1656 
Bt 660) Gh) CB cceoes cence E a wes 0. 460 
Kerosene, 
tu (ir) Gt) CB) tek oc et ee oes 1318 
BeuHseey tt) F eatencwn wend aak 0. 366 


ONE-DIMENSIONAL SPANWISE TEMPERATURE DISTRIBUTION 


The constants of integration were found by solving equa- 
tions (5) to (10) for the specific examples stated; a blade 
thermal conductivity of 15 Btu/(hr)(ft)(°F) was used. 
Blade coolant passages extending to within Me inch of the 
blade tip were considered in blades with spans ranging from 
14s to 4% inches, with water as coolant, and for effective 
gas temperatures ranging from 2000° to 5000° F in order to 


study the effect of blade length on the temperature distri- 
bution. The temperature distributions were determined, © 
by use of the calculated integration constants, from equa- 
tions (1) to (4) and are shown in figure 6. The general __ 
trend shown is a nearly constant temperature about equal 
to the coolant temperature through the rotor, a sharp tem- 
perature increase through the rim and the base of the blade, 
a nearly constant temperature (called the prevalent blade 
temperature) through the liquid-cooled part of -the -blade, 
and another sharp temperature increase to a value approach- 
ing the effective gas temperature in the uncooled part of the 
blade. From figure 6, it can readily be seen that the same 
prevalent blade temperature prevails through blades of 
various lengths cooled to within 4%. inch of the blade tip. 

In order to study the effect of coolant-passage length on 
the temperature distribution, other calculations were made, 
with water as coolant and for an effective gas temperature 
of 2000° F, for a 4X¢-inch blade span but with blade coolant 
passages of various lengths. The calculations, again deter- 
mined by use of equations (1) to (4), are plotted in figure 7; 
for short coolant passages high-temperature gradients exist 
throughout the blade, whereas for long coolant passages 
practically no temperature gradient exists throughout most 
of the blades. 
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FIGURE 6.—Effect of various blade lengths on temperature distribution for water-cooled 
gas turbine for effective gas temperatures from 2000° to 5000° F. One-dimensional spanwise 
analysis. Coolant passages extend to within Ms inch of blade tip; coolant flow, 7 pounds 
per minute per blade; thermal conductivity of blade, 15 Btu/(hr) (ft) CF). 
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FIGURE 7.—Effect of various coolant-passage lengths on temperature distribution. One- 
dimensional spanwise analysis. Blade length, 4146 Inches; coolant flow, 7 pounds per 
minute per blade; thermal conductivity of blade, 15 Btu/{hr) (ft) CF). 
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Calculations were also made for the same blades with 
4.inch coolant passages and for an effective gas temperature 
of 2000° F with water, ethylene glycol, and kerosene con- 
sidered as possible coolants... These results showing the 
effect of various coolants on temperature distribution are 
given in figure 8. For an effective gas temperature of 2000° F 
and a coolant flow of 7 pounds per minute per blade, the 
prevalent blade temperature in degrees Fahrenheit is about 
one fourth, two fifths, and one half of the effective gas 
temperature for water, ethylene glycol, and kerosene cool- 
ants, respectively. The consideration of cooling air on the 
outside of the rotor accounts for the rotor temperature being 
less than the coolant temperature. 

Finally, for the spanwise case, figure 9 shows the variation 
of coolant flow on prevalent blade temperature. The 
prevalent blade temperature decreases with increasing 
coolant flow; the rate of this decrease diminishes as the cool- 
ant flow increases. It can be seen that the cooling effective- 
ness, defined as (7, .—T7Zs)/(7,,.—T;), changes from 0.68 to 
0.91 for water as the coolant flow is increased from 2 to. 16 
pounds per minute per blade. The superiority of water is 
apparent. 
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FIGURE %.—Temperature distribution for gas turbine for various Hquid coolants. One- 
dimens{ona!l spanwise analysis. Coolant-passage length, 4 Inches; coolant passages extend 





to within 346 inch of blade tip; coolant flow, 7 pounds per minute per blade; thermal con- 


ductivity of blade, 15 Btu/(hr) (ft) CF). 
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Figure 9.—Varijation of prevalent blade temperature with coolant flow. One-dimensional 
spanwise analysis. Thermal conductivity of blade, 15 Btu/fhr) (ft) (°F). 


THREE-DIMENSIONAL TEMPERATURE DISTRIBUTION THROUGH 
as SIMPLIFIED TRAILING SECTION 

A three-dimensional temperature distribution in a turbine- 
blade trailing section approximated by a rectangular paral- 
lelepiped (fig. 2) was determined by use of equation (12) 
(approximated by use of m=1 only, as previously mentioned). 
The constants La, Qa, and N, were obtained from equations 
(16), (17), and (18), respectively; the constants A,,, and 
O,, from equations (14) and (15), respectively; and the 
constants M, n and P, a from equations (13). The distribu- 
tion was found in two planes representing the maximum and 
minimum temperatures for the Z-axis; the first plane was 
located. at the side of the cross section of the rectangular 
parallelepiped and the second plane was on the median plane 
through the rectangular parallelepiped. The results are 
shown in figure 10(a) for an effective gas temperature of 
2000° F and a thermal conductivity of 15 Btu/ (hr) (ft) (°F) 
and give the temperature distribution at various distances 
from the coolant passage. The curve labeled ‘approxima- 
tion of temperature” in figure 10(a) is a one-dimensional 
chordwise distribution through the approximated rectan- 
gular trailing section. Similar results, for thermal conduc- 
tivities of 120 and 210 Btu/(hr)(ft)(°F), are given in figures 
10(b) and 10(c), respectively. Calculation of cosMi7r/2 (see 
equation (12)) reveals the temperature variation in the twa 
planes to be about 3.6 percent of @ for a blade thermal con- 
ductivity of 15 Btu/(hr)(ft)(°F), 0.5 percent for a thermal 
conductivity of 120 Btu/(hr)(ft)(°F), and 0.3 percent for a — 
thermal conductivity of 210 Btu/(hr)(ft)(°F). Fora thermal 
conductiy ity of 15 Btu/(hr)(ft)(°F), figure 10(a) shows a 
constant temperature for the last three quarters of the blade; 


Temperature, °F 


Jemperature at cenfer of blade 
——--— femperature on surface of blade 
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Blade jes th, in. 
(a) Thermal conductivity of blade, 15 Btu/(hr) (ft) (°F). 


FIQURE 10.—Three-dimensional temperature distribution in trailing section of turbine blade. 
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Fietre 10.—Continued. Three-dimensfonal temperature distribution in trailing section 
of turbine blade. 
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FiGCRE 11.—Comparison of one-dimensional chordwise temperature distribution cecil 
rectangular trailing section with three-dimensional temperature distribution through rec- 
tangular parallelepiped. Thermal conductivity of blade, 15 Btu/(hr) (ft) CF). 





Distance from 
coolant passage that is, conduction to the rim affects about the first quarter 
ffn) of the blade length. Figure 10(a) also shows that the level 
of the temperature in the principal portion of the blade rises 
rapidly as the distance from the coolant passage is increased. 
Near the coolant passage, the distribution is in good agree- 
ment with the one-dimensional spanwise distribution pre- 
sented in figures 6 and 8. As the thermal conductivity is 
increased, as shown in figures 10(b) and 10(c), the part of 
the blade affected by rim conduction increases; for a thermal 
conductivity of 210 Btu/(hr)(ft)(°F), about two thirds of the 
blade length shows this effect. 


ONE-DIMENSIONAL CHORDWISE TEMPERATURE DISTRIBUTION 
THROUGH SIMPLIFIED SHAPES 


Temperature, °F 


Because the three-dimensional distribution resulted in a 
constant spanwise blade temperature in the region of the 
blade beyond the influence of rim cooling and because the 
critical blade point, as determined from a temperature-stress 
relation, may likewise be beyond the region of rim-cooling 
influence, one-dimensional chordwise temperature distribu- 
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FieCRE 10.—Concluded. Three-dimensional temperature distribution In trailing section ; ‘ . 
of turbine blade. tions were obtained. Inorderto compare the three-dimen- 
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FIGURE 12,—Qne-dimensional chordwise analysis of effect of shape on prevalent blade tems 
perature, Thermal conductivity of blade, 15 Btu/(hr) (ft} (°F). 


sional distribution with a one-dimensional chordwise 
distribution, a rectangular trailing section was first consid- 
ered. Figure 11 shows this comparison for an effective gas 
temperature of 2000° F and a blade thermal conductivity of 


15 Btu/(hr)(ft)(°F); the one-dimensional distribution was 


determined by use of equation (19). From figure 11 it can 
be seen that the trailing-edge temperature increases from 
1500° to 1850° F as the distance from the blade root is 
increased from 0.25 to 1.0 inch for the three-dimensional case, 
even for low-conductivity material. On the other hand, in 
the region of the blade beyond the influence of rim cooling 
(when the distance from the blade root is 1 in. or more), the 
one-dimensional and three-dimensional results approach 
identity; for such a region, a three-dimensional solution is 
unnecessary when low-conductivity materials are considered. 

In order to more nearly approximate a blade trailing sec- 
tion, a trapezoidal approximation was considered. For 
comparative purposes the trapezoidal and rectangular sec- 


lions wére constructed to have equal lengths and areas. In 


general, however, the trapezoid is constructed so that tbe 


thickness at the coolant passages equals one-half the actual 


blade surface exposed to the coolant. This dimension, in 
turn, fixes the length of the trapezoid. 


The. temperature distribution through a trapezoidal sec- 
tion was obtained by use of equation (20). Such a distri- 
bution, for an effective gas temperature of 2000° F and a 
blade thermal conductivity of 15 Btu/(hr)(it)(°F), is 
compared with the distribution through a rectangular sec- 
tion in figure 12. The temperature distribution for the 
trapezoidal section has a slightly steeper slope than that for 


the rectangular section at distances remote from the coolant _ 


passages ‘and the temperatures at the trailing edge and the 
coolant-passage wall are lower. Part of this lower tempera- 
ture for the trapezoidal section ie due to the additional thick- 
ness of the cooling surface. 

The effect of varying thermal conductivity from 15 to 210 
Btu/(hr)(ft)(°F) on the temperature distributions along the 
center line of rectangular and trapezoidal sections is shown 
in figures 13(a) and 13(b), respectively. In each case, as 
thermal conductivity is increased, the temperature distri- 
butions. flatten and approach linearity, decreasing in the 
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(a) Rectangular trailing section. 
FIGURE 13.—-One dimensional] chordwise analysis of prevalent blade temperature for trafling 
7 section of blade. 
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FIGURE 13.—Concluded. One-dimensional chordwise analysis of prevalent blade tempera: ee ae a 
ture for trailing section of blade. 90 o; 
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trailing section and increasing near the coolant passages. 


Figure 14 shows blade trailing-edge temperatures for various : Sn een nie vee a eae en ne n 
e >» > > IGUERE = railing e ratūres o om on mensio c wise analys 
thermal conductivities obtained by use of the rectangular i scene 


120 [80 £00 240 
Thermal conductivity, Btu/(hr) (77) (°F) 


Coolant temperature, 200° F. 
and trapezoidal approximations. Cooling is substantially 


greater for the more representative trapezoidal section. For 
currently used high-temperature alloys with thermal 
conductivities in the neighborhood of 15 Btu/(hr)(ft)(°F), 


; a A — /.188". 
figure 14 shows a trailing-edge temperature difference for 


the rectangular and trapezoidal sections of only about 60°F Gas-to-blade heat- ae P 

for an effective gas temperature of 2000° F. ee fea cael 1 
Chordwise temperatures were also calculated, by use of pee (Efuf(hr) (sq ft} (°F)) 

equations (22) to (27), for concentric circle annuli and 370~ i 220 “SBE 

parallel-plate approximations. These temperatures are not 380- 

plotted in this report; they were only used as initial approxi- 390- + 

mations for the two-dimensional numerical calculations made 320.: 


for an actual blade shape. 300-- 
TWO-DIMENSIONAL CHORDWISE TEMPERATURE DISTRIBUTION 

Two-dimensional temperature distributions were deter- 
mined for an actual blade shape by application of the relaxa- 
tion method. Separate calculations were made for the cases 
where an average gas-to-blade heat-transfer coefficient and 
where a typical variation in the gas-to-blade heat-transfer 
coefficient as shown in figure 15 were considered. Initial 


260-+ 6-250 


FIGURE 15.—Cross section of turbine blade showing typical variation of beat-transfer coeffi- 
cient around blade perimeter. 
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trial solutions were determined by use of equation (20) for 
the trailing section, equations (22) to (24) for the leading 
section, and equations (26) and (27) for the thin-wall sections 
of the blade near the coolant passages. A comparison of 
the temperature distributions for assumed variable and 
constant gas-to-blade heat-transfer coefficients, for an 
effective gas temperature of 2000° F, an average coolant 
temperature of 200° F, and a blade thermal conductivity of 
15 Btu/(hr)(ft}(°F) is shown in figure 16. The blade 
temperatures obtained are nearly the same for both cases 
except at the leading and trailing edges. Use of the average 
coefficient gives a conservative estimate of the trailing-edge 
temperature and a temperature that is somewhat too low 
near the leading edge. B 

A calculation has also been made for the average gas-to- 
blade coefficient and for a blade with a thermal conductivity of 
100 Btu/(hr)(ft)(°F). A comparison of this solution with the 
similar one for a thermal conductivity of 15 Btu/(hr) (ft) (CF) 
is shown in figure 17. The high-conductivity blade has 
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FIGURE 16.—Effect of variation of gas-to-blade heat-transfer coefficient on temperature 
distribution through crosa section of water-cooled turbine blade, Thermal conductivity 
of blade, 15 Btu/(hr) (ft) (°F); effective gas temperature, 2000° F: average water tempera- 
ture, 200° F. 
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FIGURE 17,—Effect of variation of thermal conductivity oi temperature distribution through 
water-cooled turbine blade. Gas-to-blade heat-transfer coefficient, 222 Btn/(hr}(sq fte E3; 
effective gas temperature ,2000° F :average water temperature, 200° F, 


about a 250° and a 600° F lower temperature than the low- 
conductivity blade at the leading and trailing edges, respec- 
tively. Little temperature difference is obtained in the 
center of the blade, where extremely good cooling prevails. 

The results of the various relaxation solutions are in good — 
agreement with the one-dimensional chordwise approxima- 
tions that were used to start the numerical solutions; that 
is, good. representative temperatures are obtainable by use 
of one-dimensional chordwise approximations. . 

In order to determine the effect of distance from the cool- 
ant passage on the trailing-edge temperature, trailing-edge 
temperatures were determined for the liquid-cooled blade 
with five coolant passages shown in figure 18. Various 
length trailing sections were obtained by successively reduc- 
ing the length of the trailing section shown in figure 18. The 
temperatures were determined by use of equation (20) and 
are shown for various thermal conductivities in figure 19. 
The trailing-edge temperature is reduced almost linearly as 
the length of the trailing section is decreased. The effect of 
thermal conductivity also decreases as the trailing-section 
length is decreased. 





FIGURE 18.—Cross section of water-cooled turbine blade showing location and size of five 
coolant passages used in determination of effect of varying tralling-section length on 
trailing-edge temperature, 


NONDIMENSIONAL CHARTS 


The availability of several nondimensional charts, to be 
subsequently discussed, eliminates the necessity for some 
numerical calculation. 

The prevalent blade temperature is given by equation 
(2a) as 

T» hdoT r etha hi sly 
ie holot hi, oli. 





o` (2a) 


After division by A, 2;,2, this equation may be written 
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B= T EN 
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FigcrRe 19.—Trafling-edge temperature as function of length of tralling section of blade 


for water-cooled turbine blades with five coolant passages and various thermal conduc- 
tivities. 


Subtraction of T,. from both members of this equation 
leads to 





Tee — Tep 1 
Tee—Ti 1A) 
Tz e—Te.p . a ° 
The plot of “or. against À is shown in figure 20. 
z.¢ +i 


For any given blade, ` can be evaluated and the value of 
Tre T B, p 
Tre- T. L 
gas and coolant temperatures, a single simple algebraic oper- 
ation results in the desired value of the prevalent blade 

temperature. 


can be obtained. Finally, for the desired effective 


Another nondimensional chart, which gives the one- 
dimensional chordwise temperature distribution through a 
rectangular section, is also presented in figure 21. The 
temperature distribution through a rectangular section is 
given by equation (19), 


i (T,.e—T1) cosh gy’ 
Le—-La= z (1 9) 
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FicgURg 20.—Nondimensional chart for use in determination of prevalent blade temperature. 


For a given turbine and set of turbine operating conditions, 
values can be determined for all the quantities in this equa- 
tion. A semilog plot with 97’ as abscissa and h./yks as a 
parameter results in a nondimensional cooling ratio 





Tae— Tsg 
(T, e— Th cosh gy’ 
The addition of a second quadrant, with cosh gy’ as param- 
eter, yields values of the temperature ratio pe at any 
g,¢ +2 


point in the rectangular section. This chart is given in 
figure 21. The chart is used as follows: A vertical line is 
constructed through the calculated value of the abscissa gj’ 
and is extended to the calculated value of the parameter h,/ykz. 
From this point, a horizontal line is drawn, extending into 
the second quadrant, and intersecting several lines repre- 
senting various values of cosh gy’. Vertical lines from the 
intersection points to the abscissa in the second quadrant 


T, s— Ts 


give values of the temperature ratio p- — at various 
1, ef; I 


positions in the rectangular section. Values of the tem- 


. perature Tp are then easily obtainable. 


In order to illustrate the use of figure 21, a rectangular 
section with dimensions equal to those previously tabulated 
will be considered. For a blade thermal conductivity of 


1&2 REPORT 1066—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 




















h 


AMARA 
USB Un S 


o 
E 


(Tee ~ T)/ [Te -Tosh py’) 


JS 


ANAT 
HH 








SSE 
5 


SNY T 
NINNAS Ras 
SAAS NEE 


NN Ri 
ae SAY Sk 
es ee AME 
ae a ee VAY! 
= TNS 






























oN AN 
ENAA SAT 


ANY i 


UE jean YNN 
o 4/0 J|] I AXIS N 
a EEANN 
EHH FHA 
THAN 


4 68l 45 




















0 Of = Sera 
l ONO Pt 


FIGURE 21,—Nondimensional chart for use in determination of temperature distribution through rectangular tralling section of cooled turbine blade. 


15 Btu/(hr)(ft)(°F) and outside and inside heat-transfer 
coefficients of 222 and 2370 Btu/(hr)(sq ft)(°F), respec- 
tively, the parameters used in figure 21 are evaluated as 
follows: 





o (2h\4 2 (222) 7 E 
da F I= 15 (0. Biot | (0.055)= 
hy o 2370 ng 


pka (54.4) (15) 


The point represented by these parameters is given by the 
circle in the first quadrant of figure 21; the square on fig- 
ure 21 designates the corresponding point for a thermal con- 
ductivity of 210 Btu/(hr) (ft)(°F). A horizontal line pass- 
ing through the circle and intersecting the family of lines in 
the second quadrant of figure 21 gives as the abscissa of the 


second quadrant the values of the temperature ratio ET ia 
d] 
for various positions in the rectangular section as follows: 
T,, e—— T, 


- ~ aae rá f= 
TZT, 0.075 for gy’=0 





=0,115 for gy’=1 
= (0.285 for oy’=2 
=0.75 for gy’=3 
For T,,.=2000° F and 7,=200° F, it follows that the blade- 
temperature calculation results in the following values: 
T'3= 1865 for yy’ =0 (trailing-edge temperature) 
=1793 for gy’=1 


=1487 for gy'=2 © eee 


= 650 for gy’==3 (temperature at coolant passage) 


These results compare favorably with the calculated distri- 
bution shown on figure 12. 


CONCLUDING REMARKS 


Analyses have been presented for obtaining spanwise 
temperature distributions near a coolant passage, chordwise 
distributions (for both approximated and actual shapes) in 
regions where conduction to the rim is inappreciable, and 
three-dimensional distributions for approximated blade 
shapes. “Numerical examples based on specific blade con- 
figurations and heat-transfer coefficients available from 
unclassified sources have been presented. Although the 
analyses are exact, the numerical values used in the calcu- 
lations may not necessarily be the same as those for cooled 
turbines. The numerical examples have been presented to 
indicate the range of applicability of the various analyses 
and to present the general nature of temperature distribu- 
tions in liquid-cooled turbine blades. In the following 
paragraphs the more important results are reviewed to show 
their general guidance for design studies. 

The three-dimensional temperature distribution includes 
conduction to the rim, whereas the two-dimensional and 
one-dimensional chordwise distributions do not. Conse- 
quently, the temperatures obtained from a three-dimensional 
investigation are less than those otherwise obtained. As a 
result, it is advisable to determine the simplified solutions 
first; if the temperatures that result are not excessively high, 
a three-dimensional investigation is unnecessary; whereas if 
the resulting temperatures are large, a three-dimensional 
investigation can be made. 
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The two-dimensional solution just referred to is of necessity 
a numerical solution (because of the varying boundary con- 
dition caused by the shape of the blade). It has been shown 
by the calculations previously presented that for uniform 
outside heat-transfer coefficients, simplified one-dimensional 
distributions were in excellent agreement with the two- 
dimensional relaxation solution. When variable outside 
heat-transfer coefficients were considered, the simplified 
solutions indicated optimistic results and the two- 
dimensional relaxation solution appears essential. Increasing 
blade thermal conductivity results in raising the cooling 
surface temperature and lowering the trailing-edge tempera- 
ture and thus in a more uniform blade temperature. 


One-dimensional spanwise distributions proved to be valid 
near the coolant passages. The temperature of the cooled 
part of the blade (prevalent blade temperature) was inde- 
pendent of blade length. If the inside heat-transfer coeffi- 
cient is increased, the prevalent blade temperature decreases; 
it was shown for a particular case that the prevalent blade 
temperature for water as coolant was only about half that 
for kerosene as coolant. 


Lewis Furcut PROPULSION LABORATORY 7 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
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SYMBOLS 
The following symbols are used in this report: Subscripts: 

cA area of metal, sq ft a air 
b blade span, ft av average 
Cia. Cy constants of integration B blade 
D diameter of circle used to approximate leading | e effective 

section of blade, ft g gas 
d ratio of distance between net points 0 and 5 | 7 inside (inner with D) 

and net spacing l liquid 
A, ) Hankel functions (special kinds of Bessel e ee eee 
cae summation indices 
tH, function) n 
h heat-transfer coefficient, Btu/(sec)(sq ft)(°F) | o outside (outer with D) 

or Btu/(hr)(sq ft}(°F) p prevalent 
J : ifs rotor 
A Bessel functions p Ae pee 
j chordwise distance from blade trailing edge to | rim rim 

coolant passage, ft T blade tip 
K } Í denote sections in spanwise investigation 

kesin y, 2 when used with x; land 2 denote trapezoidal 

k thermal conductivity, Btu/(sec)(ft)(°F) or | 8 thicknesses at trailing edge and coolant 

Btu/ (hr) (ft) (°F) 4 passage, respectively, when used with r 

perimeter, ft 4 denote end near blade tip and end near rotor 

Vf, VN, 0, P points on figure 3 B hub of various sections 
Q) heat flow, Btu/sec Superscript: 
r radial distance from hub of turbine, ft prime linear dimension inereased by 7/2 
T temperature, °F Functions: 
r distance from blade tip to blade element, ft A Gas hola \3 
y distance from trailing edge to blade element, ft kadn, 
Z number of blades B ho waned Meals 2li2\4 
3 distance from median plane of section to blade Ai g 

element. ft D (Piret hasha olos ete abs a A 
F MS Sot Fe) Ftd (hte F 2) ee | fee ee i 
t 2K| y! ya 2 Y) $ E l TA 

ney i: 4TTs ar ha Tat (hi sli da Tr] 

g T a Ts | eae Erim Ártm Arim | 

loto F\2 
A hialin G (£) 
g distance between parallel plates, ft anny, ai th, alie, az |} 
T one of trailing section, ft H | Eseta t sal der Pree eee de] 

I Feral T. A ihi. Paalse T j 

s (Fe BT. kr Ar 


7 (a) 
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K, L, M, N, integration constants (See equations (13) to 
0, P, Q (18).) 
R keM,,., sink Mya?’ +h, cosh Mun. ng’ 
4h, py 
s, HTT T 


Tx 2(T,, 





sin Q,7’ 1 
a E ) (jum) 
“20,7” 


APPENDIX B 
ONE-DIMENSIONAL SPANWISE TEMPERATURE DISTRIBUTION 


The heat balance for the uneooled_section of the blade 
(fig. 1, section 1) is given as follows: 


Heat entering by radial conduction = E z = 
l 





Heat entering sides by convection=/,/,(T, e — Te pdz 
Heat leaving by radial conduction 


=—ksAn 1 (Ts, Ged dz,) 


The heat-balance equation is 














—hpAp oa “+h LT, e— Tp, Ddr, 
=— be An A (1, +S Edn) 

or 

CT 

TIa AT, =A Tp (B1) 
where 

holo \4 
a=( 5) 
A solution is 
Te i= Tp — C cosh A{z,+ Ci (1) 


where C, and C; are integration constants. 
For the cooled section of the blade (fig. 1, section 2). the 
heat balance is given as follows: 
Heat entering by radial conduction=— ksAs,?2 tias 
2 
Heat entering sides by conyection=A,/,(T, ..—Ts..)d22 


Heat leaving by radial conduction= 


—kzAg, 2dr, (Ts Had ; dza) 
Heat leaving by radial E Ai, (TBa2— Ti) dz 





The heat-balance equation is 


—hkzApg 2 are : 





th, L(T,.—Ts, a) dx» 


=—hpAp dz, “(Ts ah ae 2 dita Hha, oli, {TB — Ti) dz 


or 
ie _ BT, =D? (B2) 
where 
pò “eet has oth, alt, 2 + 
kyApe 
and 


D=(*#74 alo Ts, et hi, ala, ra) 
ksp, 2 
A solution is 
Ta, 2= Tg, p— Cye8*2— Ce Ba (2) 


where C, and C, are integration constants and 


T»»=(5 J- 


For the rim section (assumed constant area) of the rotor 
(fig. 1, section 3), curvature was neglected because the rim 
thickness was small in comparison with the rim radius. In 
addition, average values of hala and rs were used. The 
heat-balance equation reduced to 


he lo T, eth ali, gl’; 
 AloF ha, oli. 





ERE T w= F? (B3) 
where 
E:= ÅT Fs av aha + (ha, ala, aav 
FA 
and 
pr tl s aoha Ta t (As, als, sles ly 
E rim< irim 
A solution is | 
lac G+ Cet O7,e— Bs (3) 


where C, and C, are integration constants and where 


2 
o- 

For the rotor section (assumed constant strength) (fig. 1, 
section 4), average values are used for hy alra and for r, and 
an approximate solution (as in section 3) is obtained. Be- 
cause the rotor is liquid-cooled, such a solution is adequate. 

The resulting differential equation for this section is 
d’Tpr 








JA 7 —HT,=—P (B4) 
where . 
3 TT 4 av ha+ (he, als, )av 
A Sy 
and 
pre tt ahelat (hs ali dadi 5 
rAr 
A solution 1s 
=J+C;, cosh H, (4) 


where C; is an integration constant and where 


(4) ia 


The boundary condition 


has already been applied. 
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APPENDIX C 


THREE-DIMENSIONAL TEMPERATURE DISTRIBUTION THROUGH SIMPLIFIED TRAILING SECTION 


A three-dimensional temperature distribution through a 
rectangular parallelepiped, an approximation for the trailing 
section of a liquid-cooled turbine blade, is given as follows 
(fig. 2): 





Heat entering element from top= —kgdy dz Sia 
; ; OTs 
Heat entering element from right end=—kgdrdz Sy 
OTs 





Heat entering element from front=—kgdzdy ae 


Heat leaving element at bottom 


Heat leaving element at left end 
=—ksdz dz ($2457? du) 


Heat leaving element at rear 


























=—kgdiedy oe T a dz ) 
The heat-balance equation is 
—kedydz 1? —kydrdz 28 kpde dy 22 
—kadydz (a Ned k sdzda (S84 r sed y)- 
ksdzdy (5+ ar 2) 
or r 
ols B O Tp Ò Tp 
ar oa nba B. (C1) 


Simplifieation in applying boundary conditions (as fully 
described on p. 4) and use of the effective gas temperature 
T,,¢ instead of the blade temperature Ts as the reference 
temperature, that is, use of the substitution 


6=T,.—T 
change equation (C1) into 
076 . 076 , 0768 
a ayt det 01) 


The origin of the coordinates chosen is shown in figure 2; 
the plane z=0 is the median plane of the rectangular para- 
lelepiped and from considerations of symmetry no heat flows 
across this median plane. The boundary conditions to be 
upplied are 

Og 


5,79 when 2’=0 (C2) | 
O06 ME . 
a” when y’=0 (C3) 
of =0 when z=0 (C4) 
6=T, -—Ts, when x’=0’ (C5) 


272483—3i——18 


OT, , 07T; 
=—ksgdydz al a 24 "a dz) ‘where K, L, M, N, O, P, and Q are constants. 


kag =h: (T,,.—T,—6) when y’=j’ (C6) 
and 
òg 
kg— =— h, Jhen z=1/2 (C7) 


Oz 


A solution of equation (11) satisfyi ing the three boundary 
conditions (C2), (C3), and (C4) is E 


é= D 2 Ku, 008 Ly x’ cosh Mm, ay’ cos Naz + 


mæl nel 


On,. cosh P, a t’ cos Qa y’ cos Nx Zz) (12) 


Relations | 
among these constants are 


Mon. nN P+ L,? 
Pa: NA+ Qa. 


The boundary condition expressed by equation (C5) leads 
to a determination of L. Substitution of equation (12) into 
equation (C5) gives 


njn=1,2,3, 3 (13) 


Ts .e— Tp, r= > 5 (E m,» cos L, 6’ cosh My, ,y’ cos Naz + 


mæl n=l 


On,. cosh Pm, a b’ cos Q, y’ cos Naz) (C8) 
It can be seen that the first member in.the double summation 
in equation (12) vanishes when cos Z,x’=0; hence, if 


ton(n-§ 


from equation (C8) it follows that: 


T,e— T= = >= On, cosh Py , b’ cos Q, y’ cos Naz (C9) 
This is a Fourier development along y’ and 2. 

The development of equation (C6) leads to a determination 
of Qa. Equation (C6) becomes, with the aid of equation (12) 
and its derivative, 


31 D kalE n, x (cos Lat’) My, sinh My, j’ cos Naz — 


mæl n=] 


O n,a (cosh Pum,» £’) Qa sin Q, 7’ cos Naz] 
=> > ħhi(Tre—Tiı—EK n,a cos Ly x’ cosh My, , J’ eos Nat — 
mæl n=l ge 


Onn cosh Pm, 22x’ cos.Q, j’ cos Nz) (C10) 


In order to simplify equation (C10) at the boundary when 
y’ =J and to solve for the constant Q,, the second summation 
in the left member is equated to the last summation in the 
right member term by term. Then 


h(Te.—T)= = Ke, „(keMn n sinh Ma aj’ 4 
Rm | tm 


h; cosh Mya, aJ’) cos L,x’ cos Nae (C11) 
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and 


> 32k3On, „(cosh Pia. atha sin Qo F cos Naz 


mml nml 


SSS Zh On,n cosh Pa, a2’ cos Qj’ cos Naz (C12) 


mmt ml 


The sums are equated term by term. 
determined as any solution of 


Therefore, Q, is 





tah Qn) =p =g a7) 


The development of equation (C7) leads to a determina- 
‘tion of N,,. 


> Sk, | Kn, „ cos Laz cosh Mn, ny (Nn sin Nin 5 )+ 


m=] n=! 


On,» cosh Pn, a2’ cos Quy’ (Nn sin Nn 5) -5 Sh 


mal nal 


ho (En. n cos Lpr’ odi Mm, ny’ cos Nm zt 


Om, n cosh Pm a’ cos Quy’ cos Nm 5 (C18) 
or 
hot 
tan Vn 5= = EWN zta a 


The possibility of determining values of K,,, and On,» to 
satisfy equations (C9) and (C11) has been established in 
textbooks on Fourier and other harmonic series (reference 7, 
pp. 118-121), and it is only necessary that the values 
determined define a convergent series. 

Values of K n, a are determined by integrating equation (C11) 
between the limits z’=0 to «’=06’ and z=0 to z=r/2 and 
substituting the values previously determined for L, and 
Nn» in equations (16) and (18). The integration is accom- 
plished in two steps using the functions cos Lx’ dz’ and 


cos N,z dz as multipliers, where s and v are integers. For 
abbreviation, 


R=ksMn n sith Mn aj’ +h, cosh Mm nd’ (C14) 
Then 
gimp! 
DRK», » COS Naz | | m cos"Z,x7’ dx’ 
=h(T, .—T;) TA se cos Lag’ dz’ (C15) 


because all the terms 
bf 
f cos L,z’ cos L,z’ dx’ 
0 


vanish if sn. Integration of equation (C15) and sub- 
stitution of limits lead to the result 


¥ 





SRK, Ne [in pa 


maj 
hf 
=h, | HeT ) ain L ar | 
0 


sın iat, | 


(C16) 


or 
>) RE,» cos N,,2 C 25,0" 


sin fab’ 


ee (es A a (C17) 

When equation (16) is used for La, there results 
a p E T Et ST ec 
2 RE n.n cos N,2== hi (2. —T) = =§, (C18) 


When equation (C18) is integrated in terms of z, using 
the multiplier cos M, z ds, it follows that 


T 
7=— 


RE mx cos? Nazdz=8, | “cos Nazdz (C19) 
7 s=0 : 


gad 


because all the terms 


T 


3 
f cos Naz cos N,z dz 
0 F 


vanish if mv. Upon integration, 


SN (Ne TO Sa Eon (C20) | 











or 
’ T . T 
enei) NE om 
N2 2N m 7 N m 
Therefore, 
i 2S, sin Vn 5 Lo E 
Naz sin 2Nn 5 R (C22) 
1+ $s. S A 
ms i 
2N m 3 








sin Nm > 
On, 2=2T, — 
N m 3 
cosh P; nO 
where 
T,=2 Cr e™— — Ty Ja T GET C24 
28,7 
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Finally, substitution of equations (C14) and (C18) in equation (C22) and substitution of equation (C24) in equation (C23) 
lead to the following values of-the coefficients Km.» and On, x: 























‘ T 
E, o Nea Sh Cat 7 ) Oo ao 
= (2n—1) Na | ee NT r ke My, sinh Maing’ +h; cosh Mua, x3 
2 i+ — 
2N y 5 
and 
: s sin Na- 
O _ Sil Qx) _ 4(T, «—Ts,;) 2 I 2 1 P $ (15) 
m, a Eg. in 9 a 7 f ° 
Qaj (14e) N.I dn 2Na cosh P m.ab 
2Qa) 2 f; ie 2 
Nn x - 
"9 
APPENDIX D 
ONE-DIMENSIONAL CHORDWISE TEMPERATURE DISTRIBUTION THROUGH SIMPLIFIED SHAPES 
Rectangular trailing section of blade removed from in- ; aana G AA i k A 
fluence of rim cooling. The heat balance for an element of Q= Eso 
the rectangular section is given as follows (fig. 4(a)): a g (D5) 
: : : de j 
Heat entering by conduction (right end)=kar Fy when 7 is not too large. Calculations revealed that the _ 
Heat entering by convection (sides)=2h,0dy difference between C, and 7/2 varied from 0.4 to 0.5 percent 
Heat leaving by conduction (left end) for the values of o considered and had no appreciable effect 
do a P 0.005 foot or more inside the trailing edge. 
=kzr( (i dy) 
i H From equations (D2), (D3), and (D5) it follows that 
The heat-balance equation Is kar Ce sinh (ej)=hır[T, —Tı— C cosh loj’) 
ker + 2o 6dy=kar tase v) and therefore 
or É (T,e—Ti) ` 
d*9_ sy ™—6=———_—_—_— (D6) 
dy * . oe : o sinh eit cosh (yj’) 
where 
=(7 i Substitution of equations (D5) and (D6) in equation (D2) 
P ker leads to the final equation 
A solution is | 
6= C's cosh 9 (y+ Co) (D2) z “ba cosh (py’) 
where Ck and C; are integration constants. The boundary : — fanaa) +% t cosh (oj? . | (19) 


conditions to be applied in the evaluation of the integration 


constants are 
Trapezoidal trailing section of a blade removed from in- 


da -| fluence of rim cooling.—In reference 8, the temperature- 
Ker Tyan Ta -Ti)=hir(Tee—Ti—6) when y=j (D3) distribution equation along the axis of a’ trapezoidal section 
is derived. The solution obtained is expressed by the rela- 


and 
dé tion 

a D O i g= Coit) + Cni Elit) 7) 

From equations (D2) and (D4), it follows that where Ci and Ci; are integration constants, J, and iH, are 
Bessel functions, and ¢ is defined as 
kalso sinh (pC) = ha Cs cosh (eC) (1—tan y) 
2 AR| af ee 5 

and therefore = ly tog tan y (D3) 7 
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where 
eee : (D9) 
~ Ky sin p 
and 
potan= (59) (D10) 


The constants Cio and Ci in equation (D7) may be evaluated 
by application of the following boundary conditions: 





75 =0 when y’=0 (D11) 
and 
kp E =T e T,—6)when y'=} (DIJ 
or 
T= 0 when ¢=¢, (D13) 
and 
ae =hiAlT,,e— T;— 0) when [=f (D14) 
3 dt 
It follows from equation (D8) that when y’=0, 
7{1—tan yyl} 
peek | Hoey) akg (D15) 
and when y/=7’, 
= 7 ri(1 — tan W) E 
a 7t- pr orem fa 7 (D16) 
Differentiation of equation (D8) gives 
dot do | 
, de 2K? dy’ (D17) 
hence, 
dg dé y 
je when aw 
From the properties of the Bessel functions, ` 
dd, 
em Si) 
and 
dH, 
pe =— Ey) 
Differentation of équation (D7) therefore gives 
do 
ade e ~C t) H Oa Hit - (D18) 
and it follows that 
Ciotdiltti 
C = Ett) eee 
"= E) — 
when 
dé 


From equation (D12), with the use of equations (D7) and 
(D17), it is found that 


ke ba : 
[— Cpt ditta) + Cui Lia ta] 


E 
=h.[T,.—T;— Crd a(t ta) — CutHeli ta)] 


Ta 





(D20) 


The values of the integration constants Cio and C,, are 
now found by solving simultaneously equations (D19) and 
(D20). Insertion of these values in equation (D7) gives the 
desired equation for the temperature distribution through 
a, trapezoidal section 





5 SHEE Tae T) EGY IE) FEA LGD) 
= 
BAGAGEN it r 
where 


r=H Gtit tiS hiri Hita 


For a wedge, r:£0. Then, from equation (D15) it is seen — 
that ¢,;=0, and asa consequence, Ailt) = œ and id (it) = 0. 
Equation (20) then simplifies to 


hifa 
OK%, tee 


it F(t $2)—tSi(t Fs) 





DiI) 


i= - (21) 


Concentric circle annulus approximation for blade leading 
edge.—The derivations of the equations giving blade teni- 


peratures at the inner and outer edges of the annular region 
are given as follows (fig. 4(c)): 


Heat entering from hot gas= on? —* z hol Ts, «— Ts) (D21) — 


Feat flowing over circle hi ae 


_ 2 dT, 
Date 5 "2(2) (D22) 
Heat. leaving to coolant=2-x D, hi(Ts:—T;) (D23) 
A solution of equation (D22) is 
_{_@ ( D 
Ta=( z) lo. (Cu 5) (D24) 


where Ci, is an integration constant. Application of the 


boundary conditions 


Ts=Ts,, when D= D, 
and 
Te= Ta : when D=D, 


to equation (D24) results in the elimination of Cy, and an 
evaluation of Q, namely, 


(2 = 
Ta, — Ts =( 375p ) log. D 
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or 
2rk all's, o— TL; B, ) 


D, (D25) 
OL. D: 


Q= 


From simultaneous solutions of equations (D21), (D23), and 
(D25) it is found that A 


Tule a 


(25) 
B 
L+ D E+tDT) 
log: 7, l 


From a simultaneous solution of equation (D24) and equation 
(D24) with the inner boundary condition applied and use of 
equation (D25), it is found that 





Tg=Tait D, Der, o— Tp, i) 
log. 


(24) 


From equations (D23) and (D25) and from equations (D21) 


and (D25) there are obtained 


9). 
To= Tit —— “2 (Tne Ta.) (22) 
D; hi SED D, 
and 
Tris Ta~ aks (Tz, o— Tp 1) (23) 
Doho log: 
i 


Equations (22) to (25) express the blade temperatures in 
terms of known quantities. 

Section of blade approximated by parallel nites —The 
derivation of the equations expressing the blade temperatures 
at the inner and outer edges of a blade section approximated 
by a region between parallel plates is given as follows 
(fig. 4(d)): 

Q=heat entering from the hot gas=A,(7,.— Ts.) 
Q=heat leaving to the coolant=/,(T7'3 t— T) 


@=heat flowing through the section=~? (Ts,o— Tg,:1) 


where o is the distance between the plates. Equating these 
heat flows results in the following system of equations: 


#207, Tp, D= ho (T;,e—T3,.)=Ai(T3,1—Ti) (D286) 


A simultaneous solution of equations (D26) gives 


niea = (he Teeth. T;) 
a j (26) 
LHir Che + hi) 
and 
Teta ae Pe th:Ty) i 
Tz. = (27) 
ih + hy) 
From equations (26) and (27) it iio that 
Tae Tam PE - (D27) 
a j, He +A) 
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